FABRICATION AND 
CHARACTERIZATION OF Pd/n-Si 
SCHOTTKY CONTACTS 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 

Master of Technology 


by 

Pragya Tripathi 


to the 

Materials Science Programme 

Indian Institute of Technology, Kanpur 
August, 1997 



Tfitsi 

fatfrf 

tlPraff 

lijsiiotl* 


-9 Str W'Jl 


CENTKol i-itsKA(W 

i. 'f KAMiHJft 


A iV3viiO 


^ p- M-TRI 





\ 



CERTIFICATE 


This is to certify that the work contained in the thesis entitled 
Fabrication and Characterization of Pd/n-Si Schottky Contacts 
by Pragya Tripathi has been carried out under my supervision and 
has not been submitted elsewhere for a degree. 



^itewdniL Kjuwflr 

August, 1997 Jitendra Kumar 

Professor 

Materials Science Programme 
Indian Institute of Technology 
Kanpur - 208016 



Acknowledgements 


I wish to express my deep sense of gratitude to Professor Jitendra Kumar 
who played the key role in bringing about the success of this work. The 
successful completion has been possible only due to his lively guidance, 
intellectual support, meticulous observations and critical analysis. 

I wish to thank Prof. K.N. Rai for allowing me to use the cryostat 
alongwith the close cycle helium refrigerator system. 

I ’m grateful to Dr. Subhash Chand for helping me in getting acquainted 
with various instruments and gtriding me in understanding the theoretical 
aspects 

I’m thankful to Ruchi, Sanjil and Renu for their invaluable help . 

Thanks are due also to my lab colleagues Madhukar, Tayyab, Lakhmani, 
D.K. and Prasad. I sincerely appreciate the cooperation given by my class- 
mates and friends throughout my M.Tcch. programme. I acknowledge the 
immense help received from all the officers and staff members of Materials 
Science Programme and Advanced Centre for Materials Science. 

Lastly, I would like to thank my family members for being a constant 
source of encouragement throughout this work. 


Pragy^Plripathi 


« t* 


fit 



Abstract 


Pd/n-Si Schottky contacts have been fabricated by depositing a 
palladium film of ~ 80 nm thickness in vacuum ^ 10~^ iii % using an elec- 
tron beam evaporation source onto a cleaned n-type <111> silicon wafer 
and subsequent vacuum annealing at 200‘^C“500®C for different lengths of 
time. Their current- voltage (I-V) characteristics have been measured in the 
tomperatme range 80K-300K and analysed on the basis of the thermionic 
emission-diffusion(TED) mechanism and parameters, viz., barrier height , 
ideality factor , saturation current etc. are determined. The zero bias 
barrier height is found to decrease while ideality factor(77) increases 
with decrease in temperatiue. Further , the changes in and r] are quite 
significant at low temperatures. This abnormal behaviour has been ex- 
plained by invoking the concept of barrier inhomogeneities and describing 
them by double Gaussian distribution of barrier heights with parameters as 
: zero bias mean barrier height {(j)ho) 0.93V and 0.63V, standard deviation 
{(To) 11.2 and 67 6mV in the temperature range 165-300K and 80-165K, 
respectively. 

It is shown that thermal annealing at 400*^0 , 450*^0 and 500°C of 
Pd/n-Si invariably leads to the formation of silicide of composition PdaSi 
possessing hexagonal structure (F 02 P typo ) with i)aramcters a=13.055A , 
c=27.49A. Also the growth occurs preferentially in certain directions such 
that (11.4) or (00.1) planes of hexagonal PdaSi lie parallel to the underly- 
ing (111) Si substrate in the beginning, but, shows polycrystallinity with 
time at all the temperatures. The reflectance is shown to usually decrease 
with the growth of PdaSi. In the intermediate stages a slight increase in 
reflectance observed in particular cases is attributed to the disposition of 
PdaSi preferentially on the underlying (111) Si substrate , As the time 
lapses , PdaSi layer exhibits polycrystallinity at all the temperatures and 



corresponding reflectance decreases - taking a constant value eventually. 

Finally, simulation studies have demonstrated the dominance of 
TED mechanism for current transport through the Pd/n-Si Schottky con- 
tacts. Also, negative temperature coefficients of mean {(l>bo) and variance 
{al) of the Gaussian distribution fimction (representing the barriei inho- 
mogeneities at the Schottky contact) are shown to have opposite effects on 
the temperature dependence of apparent barrier height . WTiile the former 
causes decrease in barrier height , the later compensates the change by 
a fixed amount. The overall effect is that Gaussian distribution function 
shifts towards lower barrier height and at the same time become narrower 
with rise in temperature, i.e., effective reduction in barrier inhomogeneities 
occms. 
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Chapter 1 

The Schottky Barrier 

1.1 IntroductiorA 

Within the framework of solid state and interface physics , the solid -solid inter- 
face is traditionally studied the most . One reason is the enormous importance of 
metal-semiconductor (MS) junctions and semiconductor hetero junctions for device 
physics. The metal-semiconductor junction , also known as Schottky diode , is an 
active device capable of giving a strong non-linear response under appropriate con- 
ditions [1-2], It is found in many metal semiconductor contacts and the imderlying 
physics forms the basis for their rectifying properties.There have been three major 
periods of industrial importance for the metal-semiconductor contacts - in early 
twentith century they were used as radio detectors , then as radar detectors during 
the second world war and later on as clamping diodes, microwave diodes and gates 
of microwave transistors [3]. 

The first metal-semiconductor device was made, during the early days of teleg- 
raphy , by pressing a sharp metallic wire to an exposed semiconductor surface,for 
use as radio- wave detector [3, “S J. But, they proved to be highly unreliable in 
their characteristics and were subsequently replaced by reproducible planar metal- 
semiconductor interface , produced by depositing films under high vacuum , with 
an uniform contact potential and current distribution over the junction area. 

1.2 The formation of a Schottky barrier 

An ideal surface essentially is one having atomically sharp profile with no in- 
termixing of the components on the either side , e.g., metal and semiconductor in 


1 



The SchoUky Barrier 


2 


case of a MS contact. Its electronic properties are described primarily by the work 
function (prn of metal , electron affinity Xs of semiconductor, position of their Fermi 
levels and band structures. When the two are brought into contact , equilibrium re- 
quires that the chemical potential on both the sides of interface should be the same. 

This means that the Fermi levels in the metal and the semiconductor be aligned. The 
matching of the Fermi levels requires charge to flow from one side to the other In 
the process, free electron concentration in the semiconductor near the boundary 
region decreases and, in turn, causes corresponding shift of the Fermi level. Since, 
the seperation between the conduction band edge and Fermi level E/’ increases 
with decrease in electron concentration and E/- becomes the same throughout at 
equilibrium, the conduction band edge Ec bends upward (Fig. l.l(i)). If the band 
gap of the semiconductor is not changed, the valence band edge will also move up 
parallel to the conduction band edge Ec- Also, the electron affinity (xs) of the 
semiconductor is assumed to remain unchanged. Hence, the vacuum level should 
similarly bend upwards to maintain the value of Xs- Thus, the vacuum level remains 
continuous across the transition region. Depending upon the work fimction of the 
metal and electron affinity of the semiconductor different situations may arise as 
shown in Fig. 1.1 (ii) 

A dipole layer is built up at the interface due to the flow of the charge ,the 
size of which can be estimated by considering the screened coulomb potential of a 
point charge at a distance ’r’[6,7], 

4>{r) = {q/r)exp{-r/rTF) (1-1) 

where q is the electronic charge and ttf = 0.5(n/ao)~^/® is the Thomas Fermi 
screening length; a,, is the Bohr radius and n is the electron density. With n~ 
10^^cm~^(i.e., for metals), rri? is ~ 0.55 A and with n~ 10^®cm“^(e.g. ,for semiconductor) 
is ~ 500A . Thus a normal space charge layer (depletion or accumulation type) is 
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Figure 1.1: Schematic diagram of band bending before and after metal- semicon- 
ductor contact (a) High work function metal and n-type semiconductor ,i.e., 4>m 
> <^s,(b) low work function metal and n-type semiconductor , (c) high work func- 
tion metal and p-type semiconductor , (d) low work function metal and p-type 
semiconductor [ 6 ] 
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formed. The amount, of band bending in the depletion region is correlated with the 
positive space charge of ionized donors via Poisson’s equation 

(fV _ -p(z) 
dz^ €€o 

wheic V is the potential, p{z) is the space charge density at a distance z from the 
interface. The positive charge is balanced on the metal side by a corresponding 
excess electronic -charge extending over a subatomic distance. The maximum band 
bending eVf, is just equal to the difference between two vacuum levels(see, e.g. Fig 
l.l(i)).Thus, 


■>Vfc = (pm — <j>s (1-2) 

where V{, is the potential barrier encountered by an electron moving from the semi- 
conductor to the metal. 

The potential barrier ecp^which. has to be overcome when an electron is excited 
from thp metal into the conduction band of the semiconductor is given by 

<Pb = (pm — (1-3) 

Using equation 1.1 and substituting A', = <ps — (p^. have 

r > 


(pb — QVb + (pji ( 1 - 4 ) 

Schottky assumed the semiconductor to be uniformly doped , giving a 
constant charge density in the depletion region. The electric field strength therefore 
rises linearly with distance from the edge of the space charge layer [4]. Mott assumed 
a thin layer devoid of charges sandwiched between a uniformly doped semiconductor 
and the metal. The electric field strength within the so called depletion region is 
constant [4]. The depletion region has few mobile carriers and hence has very high 
resistance . Consequently , a large part of the externally applied voltage drops across 
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the depletion region. If a negative voltage (V f) is applied at the semiconductor , the 
depletion region width decreases and voltage across it falls from Vf, to Vb-V/.The 
electrons now see a reduced barrier and move from the semiconductor towards the 
metal causing a current flow . Since no voltage drop occurs w’thin the metal , 
remains unaff^ected and there is no electron flow from metal to semiconductor . 
Thus MS junction exhibits rectifying properties. 

The difficulties involved with this simple approach are obvious. Firstly, 
the barrier height in practise does not increase linearly with (f>m in accordance with 
> equation (I-') and (1.^). Secondly, the real work function is not a constant but 
depends on the surface orientation and gets strongly affected by types of rearrange- 
ment of stmface atoms , e.g. , relaxations and reconstructions[6] . Here the nature 
of surface and interface comes into picture. 

1.3 Surface and Interface states 

The bulk crystal structure is broken at the surface so that the coordination 
of species becomes markedly different. For the topmost surface atoms the bonding 
partners on one side are missing in total .Therefore , their wave functions have less 
overlap with that of the neighbouring atoms.The splitting and shift of the atomic 
energy levels is thus smaller at the surface than in the bulk. According to Bardeen 
broken (or dangling) bonds at the semiconductor surface give rise to localised energy 
states called surface states. These states are distributed in the band gap and are 
characterized by a neutral level (j)c, . That means when there is no band bending in 
the semiconductor, the surface states are occupied right upto making the surface 
neutral , i.e., a flat band condition (Fig. 1.2 )['^] . The states below (j)c, are donor 
like while those above behave like acceptors . Depending upon the type of surface 
states and on the position of the Fermi level , the surface state may carry a charge 
, which is screened by an opposite charge inside the semiconductor .If an example of 
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Figvire 1.2; (a) Qualitaive explanation of the origin of surface states in the tight 
binding picture . Depending on their origin , these states have acceptor or donor 
like charging character (b) Electron energy band diagram of n-type semiconductor 
with surface states (i) under flat band ,(ii) in thermal eqmllibnum with the bulk, 
and (c) in contact with metal [4,^ 
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a n-typo semiconductor is taken , we have acceptor surface states ; their energetic 
positions being fixed with respect to Ec. Flat bands upto the surface would position 
the Fermi level far above the acceptor t 3 rpe surface states; filling them completely 
by electrons and building up an uncompensated charge density.This is energetically 
unfavourable and hence band bending is there at the clean surface even in the 
absence of metal contact [6 j . 

A different approach is via the mathematical treatment of Schokley [ g ] . Here, 
a monoatomic linear chain terminated at one end is considered and the Schroedinger 
wave equation is solved for an electron under nearly free electron approximation. 
Possible solutions , with real electronic wave fimctions , are standing Bloch waves 
within the chain , matched to the exponentially decaying tadls on the terminated 
side.Hence, the electronic energy levels are slightly modified from those of infinite 
btilk. Additional surface solutions become possible if complex wave vectors are 
allowed. The calculation then yields one single electronic surface state located in 
the gap of bulk states.Electrons in these states are localized within a couple of 
Angstroms of the sxirface plane.The results can be generalised to the 2D surface[6] 

Numerous experiments on the wide variety of MS systems prepared under 
UHV conditions indicate that the deposition of metal films produce interface states 
which determine the position of the Fermi level at the contact [6] . According to 
Heine [iC] metal (Bloch) wavefunctions tail into the semiconductor in the energy 
range in which the conduction band of metal overlaps the forbidden gap of the 
semiconductor. Fig.l.3(b) represents the situation in greater detail. The break- 
down of the periodcity at the interface introduces exponentially decaying interface 
states with imaginary wave vectors, which fills the energy gap of the semiconductor 
synxmetrically with respecf to the band edges Ec and E^ . In a simple ID model 
(heir dis^rsion curve E-4 (Fig. 1.3(a)) {tidicates the range of theoretically possible 
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Figure 1.3; Origin of metal induced gap states at MS interface (a) Dispersion curve 
of a ID model of semiconductor ; the real band structure with Ec(k) and E„(k) as 
conduction and valence bands respectively. Dt,s is density of the virtual states (b) 
Tailing of metal Bloch state into semiconductor 

states known as "Virtual induced gap states "(VIGS). In the case of free surface , 
only a single energy level exists , but in the present case ,the Bloch states of metal 
conduction band have to match resulting in a broad continuum of states in the en- 
ergy gap of the semiconductor. Or else, the continuum of metal states leaks into the 
VIGS which are themselves derived from the semiconductor band structure . Which 
of these states actually exist depends on the boundary conditions at the interface. 

Fig. 1.4 shows barrier heights as predicted by the Schottky , Bardeen models 
together with experimental findings for various metal-silicon contacts. Accordingly , 
the esqjerimental result lie in between the two extreme approaches namely Schottky 
with no interface states and the Bardeen with high density of interface states. If 
the density of surface states is high , the charge exchange takes largely between 
the metal and the surface states and the space in the semiconductor nearly remains 
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Figure 1.4: Barrier heights of Si-Schottky contacts versus metal work function. For 
comparison the predictions of the Schottky (no interface states) and the Bardeen 
(high density of interface states) model are given [6]. 

unchanged. As a result , barrier height becomes independent of the metal work 
function and is given by [4 ] 


4>h = E,- (t>o (1.5) 

Here the barrier height is said to be pinned by surface states and equation(1.4) 
is known as the Bardeen limit. The assumptions of Bardeen are, of course, incorrect , 
since the surface states of the clean stuface are strongly affected by metal deposition. 
Taking into accoxmt the interface states , the Schottky barrier diode can be described 
as follows (Fig. 1.5): 

When the metal comes into close contact with the semiconductor surface , 
they’ll form chemical bonds , The distribution of intrinsic surface states of the clean 
semiconductor surface will be changed . Charge will flow from one side to other 
due to formation of bonds which is described by the formation of a dipole layer 
of atoroic dimensions. The direction of the dipoles is determined by the difference 
between the work function and electron affinity. Formation of new interface states 
is expected in all the cases. In a more extensive band scheme of the MS contact , 
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Figure 1.5: Band diagram of a metal-semiconductor junction in which the interface 
region (width ~5 A) is taken into account explicitly 

interface region could be included with its atomic dimensions. Using the metal work 
function(^m) and the electron affinity (Xs)of the semiconductor , dipole energy (6) 
can be associated to the interface layer . The total charge balance at the interface 
now includes the space charge in the semiconductor and the charges present in the 
metal and those located in the interface states. In thermal equillibrium , no net 
current flows . It is obvious that the band bending and in turn , the Schottky barrier 
height depends on the nature of the interface states.The origin and character of 
these states need to be further investigated for better understanding of the Schottky 
barrier formation. 

1.4 Conduction Mechanisms 

The various ways in which electrons can be transported across the MS jimction 
under forward bias is shown schematically in Fig.1.6 for an n-type semiconductor. 
The inverse processes occur under the reverse bias condition . The mechanisms in- 
volve: (a) emission of electrons from the semiconductor over the top of the barrier 
into the metal (b) quantum mechanical tunneling through the barrier (c) recom- 
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Figure 1.6: Transport processes in a forward-biased Schottky contact [4]. 

bination in the space charge region (d) recombination in the neutral region(hole 
injection) [4] 

1.4.1 Thermionic emission-diffusion Mechanism 

The motion of an electron crossing over the barrier from the semiconductor 
into the metal , moving through the high field depletion region, is governed by the 
drift and diffusion processes [ it]. Density of the available states, in fact controls the 
jumping (or emission ) of electron into the metal. The two processes , the emission 
over the barrier and the drift and diffusion in the depletion region are effectively 
in series and the one which offers the higher resistance determines the current. In 
Bethe’s TE theory [7'j , it is assiimed that the barrier height is much larger than 

kT and the effect of drift and diffusion in the depletion region is negligible. Hence 
, the current flow depends solely on barrier height. Schottky in his diffusion 
theory had also included the effect of electron collisions within the depletion region 
and assumed the cairier concentration at any distance to be unaffected by current 
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flow.Both the appioaches led to the similar current expressions. According to the 
synthesis of the two theories proposed by Crowell and Szc [ . i ] the net current due to 
flow of electrons from semiconductor to metal by overcoming the potential barrier 
of height at a forward bias V is given by [4,7] 

J = A,A"T\xp(=^)[exp(^) - (1.6) 

where A** is the effective Richardson constant, Ad is the diode area , T is the tem- 
perature in Kelvin , k is the Boltzmann constant and q is the electronic charge.For 
free electrons Richardson constant is 1.12 x 10® Am“^K~^. However , the barrier 
height is known to increase with forward bias. If the variation is assumed to be 
linear , we can write [7] 


hiV) = (j>bo + lV (1.7) 

where <j)bo is the barrier height at zero bias and 7 (= ^) is positive. Substituting 
the value of ^b{y) into eqn (l.o) , we get 

I = 4exp(^)[l - exp{-^)] (1.8) 

where 

I. = AiA"T^exp(^^) (1.9) 

and 

1 = 1 -. 

V 

The parameter r) is called the ideality factor and usually has a value greater 
than unity. Further , the presence of a series resistance(Rs) causes significant voltage 

P 

drop at large forward currents.As a consequence , the current equation (l.!<) gets 
modified to 
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Figure 1.7: Field and thermionic-field emission under forward bias . 


_ _ ,q{y — IRs).., —q(V — IRs) 

I = IgCxp^ )] * 

Also, the ln(I) vs V plot begins to deviate from linearity and exhibits saturation. 


1.4.2 Tunneling through the Barrier 

This mechanism prevails when the semiconductor is heavily doped such that the 
Fermi level lies above the bottom of the conduction band and depletion width is 
very small. It may be possible for electrons with energies below the top of barrier to 
penetrate by the well known quantum-mechanical tunneling. In the case of degen- 
erate (heavily doped) semiconductor at low temperatures , the current in forward 
direction arises from the timneling of electrons with energies close to Fermi energy 
in the semiconductor . This process is known as field emission (FE) . At raised 
temperatures , electrons excited to higher energies have higher tunneling probabil- 
ity as they experience thinner and lower barrier . This is known as thermionic field 
emission (TFE) [4] (Fig.1.7). If the temperature is still raised further , a point is 
reached in which virtually all the electrons have enough energy to go over the top 
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of barrier; it is pure thermionic emission with negligible tunneling.The tunneling 
current through the barrier is of the form [4, ‘'i , 



E, = E^cothi^) ( 1 . 12 ) 

and \to is the saturation tunneling current. Eqo is defined as 

= =18.5X10'^^-./— eV (1.13) 

47r V r)if*es y 

where me*‘(=nvmo)is the effective mass of the electron and (=6^6^,) is the 
permittivity of the semiconductor, mo is the electron rest mass and N<i is the donor 
concentration in m"^. The process is FE when Eoo>>kT and is TFE at Eoo ~ kT 
. TFE also causes a reduction in barrier height by [4] 

A4> = eJv,'^ (1.14) 

II 

where stands for the voltage corresponding to the band bending. Eqn (1.1 i j 
suggests that lult vs V plot should yield a straight line with slope giving Eo and 
intercept at zero bias the saturation current To- 

1.4.3 Carrier Generation and Recombination 

At zero bias there exists a thermal equillibrium , as the rate of electron-hole pair 
generation in the depletion region is balanced by the rate of recombination. With 
the application of bias, however , np product depeirts from n,^ and a net generation 
or recombination of carriers occur depending upon the polarity. In case of a Schottky 
diode made on n-type silicon ,when the substrate is forward , biased electrons are 
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injected into the depletion region from the neutral bulk semiconductor and the holes 
are swept away from the metal . These additional electron-hole pairs combine in 
the depletion region to give a forward combination current . The recombination 
current can then be described by [4] 


- 1] 

(1.15) 

qn,AiW 

With Iro = — 

2r 

(1.16) 

and n^- \l NcN^exp{ ) 

(1.17) 


Here w is the width of the depletion region , r is the effective carrier lifetime 
within the depletion region , is the intrinsic cmrier concentration , Eg is the 
energy band gap and Nc and are the effective conduction and valence band 
density of states, respectively. 

1.4.4 Minority Carrier Injection 

Under the low bias conditions , the Schottky barrier diode is a majority carrier 
device. But, at large forward bias , minority carrier current contribution increases 
and becomes more than the diffusion current. The electrons flow from the semicon- 
ductor into the metal under the forward bias and some of the holes are swept into 
the neutral region . The current I^, caused by hole injection from metal into the 
neutral region of n-type semiconductor of donor concentration is given by [ 




NdLp 


(1.18) 


where Dp is th^ hole diffusion constant , Lp = yj DpTp is the hole diffusion length 
and Tp represents the hole lifetime in the neutral region. 
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Figure 1.8; Activation energy plots for determination of barrier height [7] 

1.5 Methods of Barrier Height Measurement 


In general , the barrier height is measured either using the current-voltage 
or the capacitance-voltage characteristics^ On the basis of thermionic emission- 
diffusion (TED) mechanism , the current through the Schottky barrier diode at 
a forward bias V is given by eqn {1.^^. The plot of In (I) vs V is linear upto a 
point , for the forward bias V in excess of 3kT/q [4] . Afterwards the plot deviates 
showing saturation due to the series resistance term. In the plot of ln(I) vs V, 
the straight line portion is extrapolated to zero bias , so that the intercept at the 
ordinate is In(Is). By substituting the values of effective Richardson constant A*’‘ , 
saturation current Ij in eqn (l,b) , zero bias barrier height can be determined 
at a temperature T for a given diode area A^.The value of <f)bo is not very sesitive 
to the choice of A** since at room temperature a 100% increase in A** causes an 
increase of only O.OlSVin (f)ho [7]. 
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V (VOLTS) 

Figiire 1.9; 1/C^ versus applied voltage for determination of barrier height [7] 

Also , the activation energy (or the Arrhenius) plot i.e. In (Is/T^) vs 1/T 
plot (Fig. 1.8) can be made using saturation current data at various temperatures. 
<j)bo can be obtained from the slope of the best fit line [4]. Here the intercept at the 
ordinate is ln(A’"' A^) where Aj is the electrically active area. The main advantage 
in this approach is that the value of Aj is not at all reqrdred for the evaluation of 

^bo • 

Second method of determining the barrier height is by measurement of the 
capacitance of the Schottky diode as a function of reverse d.c. bias(VR)(Fig.l.9). 
The charges of opposite signs are induced on the metal surface and in the depletion 
region of the semiconductor when a small a.c. voltage of few mV is superimposed 
on a d.c. bias. The relationship between the capacitance (C) and reverse bias Vjj 
is given by[l] 
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^sqNd 




(1.19) 


where e* is the permittivity of the semiconductor, Nd is the dopant concentration 
, is the distance of the Fermi level from the bottom of the conduction band and 
other symbols have their usual meanings. Thus , if 1/C^ is plotted against Yr , a 
straight line with a slope of {2! Ad ^sqN^ is obtained having an intercept on the 
voltage axis Y rq which 1/C^ is zero.The value of barrier height <j)i, can be than 
easily obtained as [4] 


0b — Vro + 0n H (1-20) 

Q 

and 0,1, the depth of Fermi level below the conduction band , can be computed 
from the doping concentration can be evaluated from the slope of 1/C^ vs 

Yr plot too. 

1,6 Importance and Applications 

All the practical applications exploit the inherent high speed advantages of 
majority carrier conduction in Schottky diodes. The most important commercial 
application of Schottky diodes is in bipolar integerated circuits as clamps and to 
a lesser extent as load resistor substitutes , diode couplers and level shifters. They 
are also used as discrete microwave transistors , optical and nuclear particle de- 
tectors, etc. The Schottky diode is preferred over the pn junction diode because of 
high forward conductance per imit area , large reverse impedance , fast switch- 
ing speed about Ins , compact size ,and relatively easy fabrication technology at 
quite low temperatures. Metal-semiconductor rectifiers are relatively immune to 
the speed limitation associated with the minority carrier storage in p-n jtmction 
devices . However, the reverse current in a Schottky barrier is quite large, 
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and is disadvantageous . Also , it is a surface device and hence sensitive to surface 
contamination and perimeter effects. 

1.7 Objective of the present work 

Extensive studies have been undertaken on different aspects of Schottky bar- 
rier diodes with the view to understand their characterstics better and find newer 
applications [3,4) The most important aspect of metal-semiconductor con- 
tacts is the process which determines the flow of charges over the barrier from the 
semiconductor to the metal and vice- versa. Detailed knowledge of the conduction 
process involved is needed to extract the values of the barrier parameters. The I-V 
cheuacteristics of Schottky diodes formed by depositing various metals such as Ir, 
Pt, Cr, Er, W ,Cu, and Ti on silicon have been studied in detail [ 22-32.]. 

The fewer reports are avcdlable for Pd/Si Schottky diodes in particular and 
contain limited information on the conduction mechanisms ,e.g., Chin et. al. [33] 
studied Pd 2 Si/Si Schottky barrier diodes on n-type silicon at 300K and on p-type 
silicon at 84 K . Werner and Guttler reported some data on PdaSi Schottky 
barrier diodes while describing the temperature dependence of barrier heights for 
various metals on n-t 3 Tpe silicon. Chand and Kumar have recently studied 

PdaSi/n-Si contacts over a wide temperature range and interpreted the results on 
the basis of thermionic-emission-diffusion mechanism, assiiming the presence of a 
Gaussian distribution of barrier heights at the interface. 

The current-voltage (TV) characteristics of Schottky barriers measured at 
room temperature do not give adequate information about the conduction process 
and the nature of the barrier formed at the metal-semiconductor interface. Details 
of many possible effects that cause non-ideality in the diode I-V characteristics , 
and , in general , reduce the barrier height cannot be ascertained. The temperature 
dependence of the I-V measurements gives a better picture. Also, Schottky barriers 
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in devices operating in space and in some other applications are subjected to wide 
range of temperatures. So, information about their electrical characteristics over 
wide temperatvue range helps in understanding the actual^of these devices as well. 

An attempt has been made here to fabricate Pd/n-Si Schottky contacts 
and study their I-V characteristics in a wide temperature range (80K-300K). Since, 
fabrication involves deposition of palladium metal film and subsequent annealing , 
identification and growth of resulting phase (s) of silicide have also been explored 
using X-ray diffraction and reflectance measurements in UV-Visible range (200- 
900nm). In addition , simulation of (1) current-voltage (I-V) characteristics have 
been carried out for various mechanisms with a view to exmine their relative con- 
tributions and dominance , and, (2) apparent barrier height (BH) as a function of 
temperature . For the later case , barrier inhomogenities are believed to be present 
at the contact and described by a Gaussian distribution of barrier heights with a 
mean and standard deviation (a^). In this, uniqueness lies in considering the 
linear temperature dependence of the mean and variance (cr^). The purpose is to 
understand the effect of various parameters in detail. 



Chapter 2 
Experimental 

2.1 Fabrication of Schott ky diodes 

The Schottky diodes were fabricated upon n-type (phosphorus doped ) silicon 
wafers of (111) orientation . The njxi^ silicon wafers had a 8-9 /jm thick epitaxial 
n layer of resistivity ~2.6 ohm-cm(Na = 3.9X10^® cm~^) over the heavily doped 
n"'' region. The wafers were first degreased with soap solution and then were thor- 
oughly cleaned with organic solvents, viz, trichloroethylene , acetone and methanol 
in succession in an ultrasonic cleaner and rinsed in deionised water and dried. 

The back contact (ohmic) was made by depositing a layer of aluminium 
on the heavily doped n"^ side of the silicon wafer after usual cleaning and etching 
in dilute HF (HF:H2O=l:10) followed by washing in deionised water .A~ 200 nm 
thick aluminimn deposition was carried out in a HHV coating unit model 12A 4D 
under vacuum ~10“® mbar, at 450" C for 30 minutes after introducing in a quartz 
tube , placed in a tubular furnace and connected to a pumping system model HHV 
VS114D. The silicon wafers were cleaned again from the epitaxial n or p layer 
(i.e. front) side in trichloroetylene ,acetone and methanol and then etched in the 
dilute hydrofluoric acid (HF:H 2 O=l: 10 ) and rinsed in deionised water for removal 
of silicon dioxide layer that usually forms on the silicon wafer during exposure to 
the ambient, aluminium deposition step and annealing. The back contact of the 
al umi nium was protected from HF by laying a wax coating. Palladium film(~80 
nm thick) was deposited by employing an electron beam evaporation source in a 
HHV coating unit model 12A 4D-SC on to the cleaned epitaxial n-layer through 
circular holes of diameter 1 mm in a stainless steel mask. Subsequently, the samples 
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were taken out of the work chamber and put inside a quartz tube as before and 
annealed in vacuum of the order of~ 3-5 x 10~®mbar at 300“C, 400^C, 450°C, 500°C 
for various durations to form palladium silicide. The silicon wafers were finaly cut 
into small pieces such that each contained a palladium silicide dot and mounted 
on to a T05 header through silver epoxy. The contact to the top metal dot was 
established through a thin gold wire. 

2.2 Setup for I-V-T measurements 

An automated setup shown in figure 2.1 was used to record the I-V characterstics 
of Schottky diodes was over a wide temperature range (80K-300K). The bias was 
applied from a Keithley programable voltage source model 230 in steps of 5-lOmV 
and the current through the diode was recorded with Keithley picoaiheter model 485 
through a computer controlled system. To ensme steady state conditions, a delay 
of few seconds was observed between subsequent measuring steps. For recording 
the tempreture dependence of the I-V characterstics in the range 80K-300K a CTI- 
cryotronics close cycle helimn refrigerator model 22C equipped with a Lake Shore 
temprature controller model 330 was employed. 

2.3 Phase identification microstructure and re- 
flectance measurements 

For identification of the phase(s) formed and studying the growth kinetics 
of silicide, palladium film was deposited on a beire silicon wafer after usual clean- 
ing and etching in hydrofluoric acid exactly in manner described in section 2.1, 
but without mask. Samples were annealed at 200°C, 400°C, 450°C and 500°C for 
various durations in vacuum~ 3-5xl0”®mbar. Rich Siefert ISO-Debyeflex 2002 X- 
ray diffractometer with CuKq radiation (A=1.5418JL) was employed to record the 
diffraction patterns. The calibration of the system was carried out with-a standard 
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automated set-up for i-v measurements 
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Figure 2.1: Block diagram of automated set-up for I-V measurements of Schottky 
diode over a wide temperature range. 
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Picoammcler 
(Kcithly 485) 

Figure 2.2: The circuit diagram for I-V measurements of a Schottky diode. 


silicon sample. The observed X-ray diffraction peaks were indexed by determining 
the interplanar spacings (d) corresponding to various planes using the well known 
Bragg’s law [<f^2d sin G = X, 9 being the peak angle. The phase(s) identification 
was carried out by comparing the results with the standard data available in the 
powder diffraction files|^lj. 

The microstructme of the samples annealed at various temperatures for differ- 
ent lengths of time was examined in scanning electron microscope JEOL JSM 840A 
in the secondary electron(SE) mode. Also, reflectance measurments at near nor- 
mal incidence (by employing a 5° spacer) were carried out for the annealed samples 
in the wavelength range 200-900nm using a Hitachi UV 150-20 spectrophotometer 
equipped with a 60 <?!> integrating sphere. 



Chapter 3 

Results and Discussions 

3.1 Forward I-V characteristics 

The forward I-V characteristics of Pd^jSi Schottky diodes fabricated on n-type 
Si (111) as described in section 2.2 and annealed at 450^C for 45 min are shown in 
figme (3.1) at various temperatures in the range 80K-300K. Clearly, they are linear 
over several order of current. Further, they progressively become straight over a 
wide current range with decrease in temperature. Such an observation suggests the 
dominance of the thermionic emission -diffusion mechanism for cmrrent transport 
. A program in Basic is used to fit ^he experimental I-V data in the thermionic 
emission-diffusion current equation (1.10). It involves iteration method and takes 
Is, rj, and Rj as adjustable parameters. 

3.1.1 Zero bias barrier height , ideality factor and series 
resistance 

The zero bias barrier height (^^jo) can now be evaluated from the I^ data by 
substituting the values of the effective Richardson constant A**, diode area and 
temperature T. Thus is found at each temperature from equation (1.9). Figure 
(3.2) shows the zero bias barrier height as a function of temperature. It clearly 
depicts that <i>ho decreases with fall in temperature steadily in the begining but 
quite rapidly below 120K.The vzuiation of ideality factor and series resistance with 
the temperatxue are shown in figure (3.3) and (3.4) respectively . The ideality factor 
increases rather slowly with decrease in temperature but quite significantly at lower 
temperatures. There is a rapid increase of series resistance with fall in temperature 
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below 200K which is beleived to be due to factors responsible for increase in r) and 
/or lack of free charge carriers at lower temperatures [22.]. 

3.1.2 Activation Energy Plot 

Equation (1.9) can be rewritten as 

= (3.1) 

Therefore, the ln(Is/T^) vs 1/T ( i.e. activation energy or Arrhenius) plot should 
exhibit a straight line with slope giving the zero bias barrier height {(f/bo) and inter- 
cept at the ordinate determining the Richardson constant (A**) itself for the known 
diode area A^. This method gives a single value of barrier height. The experimental 
results , however correspond to the ln(Is/T^) vs 1/T plot , shown in figure (3.5). 
The data fits to a straight line only at high temperatures , i.e. , above 180K. The 
slope of the curve is infact varying continuously. 

3.2 Concept of Barrier Inhomogenities 

According to Chand and Kumar [ ] an abnormal decrease of barrier 

height (BH) and increase of the ideality factor (77) with the decrease in temperature 
cannot be explained adequately by incorporating interface states or interfacial oxide 
layer, tunneling , image force lowering and / or generation - recombination effects. 
A way to explain the abnormal behaviour of barrier height and ideality factor is 
to assmne the existence of barrier inhomogenities and describe them with some 
suitable distribution function. Numerous evidences have since been reported for 
the existence of barrier height inhomogenities [^#2.-^]. Inhomogenities may be due 
to variation in thickness and composition of the silicide layer, nonuniformity of the 
interfacial charges and locally defective hot regions [O]. 

Various types of distribution functions have been proposed to describe 
inhomogenities, e.g., a Gaussiab or a log fiormal. One can assume a Gaussian 
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Figure 3.1: Current - voltage characteifctics of Pd/n-Si Schottky contacts in the 
temperature range (a) 300K-180K (b) 165K-80K 
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Figure 3.2; Variation of the zero bias barrier height (^jo) with temperature ,in the 
range 80K-300K 
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Figure 3.3: Temperature dependence of Ideality factor for Pd/n-Si Schottky con- 
tacts ,in the range &PK-300K. Noticie the sharp increase in rj at low temperatures. 
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Figurp 3.4: Temperature dependence of series resistance (Rj) for Pd/n-Si Schottky 
contacts ,in the range 80K-300K. Notice the sharp increase in Rs at low tempera- 
tures. 



Figure 3.5: The ln(R/T2) iqoO/T plots of Pd/n-Si Schottky contacts with best 
fit straight line ,in the range 80K-300K. It shows non-linearity below 165K 
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Figure 3.6: The modified activation energy plots for cr^ 0.1109 V and 0.063 V 
, together with the experimental data. The straight lines 1 and 2 indicate best 
fitting of data in temperature ranges 165-300K and 80-150K. 


distribution of barrier heights with a mean value of (j)b and standard deviation a , 
having the form[ 2 - 5 ,*/r,^^ 

where preexponential term is the normalisation constant . The total current 
across the Schottky contact at a forward bias is then given by 

1{V) = r l{<j>b,V)p{cl>b)dcj>b (3.3) 

J—oo 

where l{<j)b , V) is the current for an elementary diode having barrier height (j)b 
at bias V and is given by equation (1-10) in case of thermionic emission-diffusion 
process. Substituting 1(^6 , V) and p{<j>b) from equations (1-10) and (3.2) in (3.3) 
and performing the integration from — oo to oo one gets 
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The barrier is known to depend on the electric field and hence on the applied 
bias [4] . Now , since barrier heights are represented by a Gaussian distribution, the 
entire profile should be affected by the applied bias. Hence considering the linear 
bias dependence of the distribution parameter, we have 

^b = <f>bo + l/V and a = (To + (3.5) 

where 7 and ^ are positive derivatives , and and cto are mean barrier height 
and standard deviation at zero bias, respectively. Making these substitutions and 
neglecting the term involving equation (3.4) gets modified to 

I{V) = (3.6) 

with 

l. = A^A-T^exp{-^) 

Here is the observed saturation current , (j)ap and rjap are the apparent barrier 
height and ideality factor respectively, ^op and r)ap are given by 

<l>op = ^bo ~ 2 ^ = (1 ~ 7) + (3-7) 

It is clear that <j)ap depends upon the distribution parameters (i.e. mean d^o and 
standard deviation cTo) and temperature T. Obviously, decrease in the barrier height 
is caused by mere existence of distribution and its effect becomes more and more 
pronounced at low temperatures. On the other hand, is a function of rate of 
change of both the mean (f>bo and the standard deviation cr^with the bias besides <To 
itself. 


3.3 Interpretation of the Experimental Data 


The diode parameters derived from the fitting of the experimental data in 
equation (3.6) can therefore be interpreted using equation (3.1) and (3.^) . Thus 
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thf' <f>ap VS 1/T plot should be a straight line giving <j>bo and <Tc from the intercept 
at the ordinate and the slope respectively. The I-V data however, reveal , two 
straight lines in vs 1/T plot giving two values of <f)bo {0.93V and 0.63V) and cr^ 
(0.118 and 0.067V). Such an observation is , in fact attributed to the presence of 
two Gaussian distributions of barrier height at the interface and discussed in length 
later in section 3.4 

Also combining equations (3.6) and (S."?) one obtains 

Hp) - ^ ^ (3.8) 

The nonlinearity in the activation energy plot can now be understood on the 
basis of equation (3. 8 )■ The factor In(^) — is calculated to modify the experi- 
mental data and plotted as a function of 1/T to obteiin a straight line corresponding 
to a single activation energy around (j>bc accurately. The intercept of the straight 
line at the ordinate gives the Richardson constant itself. Figure (3.6) shows the 
modified activation energy plot . Notice that the deviation is removed , but, the 
data fits with two straight lines. These correspond to the value of 4)^ 0.93 and 
0.63V and Richardson constants 0.99 x 10® and 1.69 x 10® Similarly, the 

temperature dependence of ideality factor {fjap) can be tmderstood on the basis of 
equation (3.7). It indicates that l/rjap vs 1/T plot should yeild a straight line with 
slope giving CToq^/k and intercept (1-q)- The increase of rj with fall in temperature 
arises due to the term containing ^ (i-e., the rate of change of <t with bias) . The 
plot of l/r)ap vs 1/T also contains two straight lines . Both correspond to negative 
slopes making ^ essentially negative;the values being -2.38mV and-0.14mV in the 
temperature range 80-135K and 150-300K, respectively. 

It may be mentioned that the studies on many other simultaneously fabri- 
cated diodes revealed slight variation in the values of the mean barrier height and 
standard deviation at zero bias. Such observations are primarily attributed to the 
local chemistry and conditions of the interface. 
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3.4 Evidence for a Double Distribution 

Figure (3.7) shows yet another Inl vs V charactenstics of a Schottky diode fabri- 
(•ato(i after annealing at 400“C for Ih in vacuum. These cle2irly depict increasingly 
liiK'ar behaviour over several orders of current and shift towards higher bias side 
with decrease in temperature . The variation of apparent barrier height and ideality 
factor with temperature as obtained in this case are given in figures (3.8) and (3.9) 
respectively. According to equation (3.7) , the plot of (j)ap should be a straight fine 
with intercept at the ordinate determining the zero bias mean barrier height <f>bo 
and the slope giving the zero bias standard deviation ao- However, when the (j)ap 
vs 1/T plot is made (figure (3.10)) , data fits nicely , as in the previous case , with 
two straight fines (transition occurring ~ 135K), implying thereby that (f>i,o and ao 
assume different values in the temperature range 165-300K and below 165K. The 
intercepts and slopes of these straight lines yield values of and ao as 0.93V and 
0.112V in 165-300K and 0.629V and 0.067V below 165K. (f^ap can now be estimated 
from equation (3.7) itself over entire temperature range 80-300K for each set of <j)ho 
and ao values. This is presented in figure (3.11) by continuous curves. 

Similarly , l/j^ap vs 1/T plot should also possess different characterstics 
in the two temperature ranges if contact contains two barrier heights distribution. 
This is indeed the case as data fits with two straight lines in the temperature ranges 
300-165K and below 165K ( see ,e.g. , figure (3.12)). Their slopes yield qa-o^/k and 
intercepts at the ordinates give (I-7 ) . The value of | extracted firom the slopes 
are -2.4mV and -0.14mV whereas the values of 7 obtained firom the intercepts are 
-0.75 and 0.187. 

Figure (3.13) shows the conventional activation energy ln(Is/T2) vs 1/T 
plot . The experimental data are seen to fit asymptotically with a straight line 
at higher temperatures only. The zero bias barrier height obtained firom the slope 
of this straight fine comes out to be 0.34V . Likewise the value of Richardson 
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Figure 3.8: Variation of the zero bias barrier height {4>i^ with temperature ,iii the 
range 80-300K. 
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Figure 3.9: Temperature dependence of Ideality factor for Pd/n-Si Schottky con- 
tacts ,in the range 80-300lt. Notice the sharp increase in t) at low temperatures. 
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Figure 3.10: The barrier height {<j>bo) obtedned from I-V measurements as a function 
of inverse temperature.Notice the transition occuring arotmd 150K . 



Figure 3.11: The zero-bias barrier height as a function of temperature.Full curves 1 
and 2 correspond to <j)ap calculated using equation (3.8)for two distributions having 
(j>bo (0.93V and 0.629V ) and or„ (0.112V and 0.0667V;. 
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Figure 3.12: The inverse ideality fzctoT{l / rjap) versus inverse temperatinre (1/T) 
plot. The data show linear variation in the two temperature ranges with transition 
around 200K 



Figure 3.13: Conventional activation energy ln(Is/T^) vs 1000/T plot. Notice the 
deviation from linearity below 180K. 




1000/T (K"^) 

Figure 3.14: The modified activation energy plots for cr^ 0.112 V and 0.067 V 
, together with the experimental data. The straight lines 1 and 2 indicate best 
fitting of data in temperature ranges 165-300K and 80-150K , respectively. 

constant determined from the intercept at the ordinate is 3.34x10 Am~^K“^ . Both 
these parameters are quite unreasonably low as the reported value of the barrier 
height for Pd 2 Si/nSi Schottky diodes is 0.735V [4F] and Richardson constant for n- 
type silicon is 1.12x10® [4] . Moreover, the activation energy plot shows 

nonlinearity at low temperatures. To explain these discrepancies the I-V data are 
interpreted in terms of equation (3.10) and a modified activation energy plot (figure 
(3.14)) is made. For this the term In(^) — is calculated for both the values of 
cTj, ( associated with the Gaussian distributions of barrier heights) found above and 
then plotted as a function of inverse temperatxire . The best linear fitting yields 
zero bias mean barrier heights as 0.969V and 0.680V which match exactly with 
the mean barrier heights obtained from the <j>ap vs 1/T plots earlier. The intercept 
at the ordinate is just ln(Ad A*”) and gives Richardson constant A’'* as 1.22x10® 
Am'^K"^ and 1.69x10® Am~^K~^ for the two straight lines drawn. 

All the considerations above suggest that the obsersved experimental data can 
be satifactorily exjjlained by assuming the existence of two Gaussian distributions 
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of barrier heights in the present PdaSi/nSi Schottky diodes. 

3.5 X-ray Diffraction Studies 

Fig. (3. 15), (3.16), (3.17), (3. 18), and (3.19) show X-ray difl&ractograms of grow- 
ing silicide at different stages of annealing at 200°C , 400°C ,450‘'C and 500^0 . 
These are recorded with CuKa radiation(wavelength A = 1.5418 A ) at the same 
sensitivity and time constant. The Bragg angles (20) , d- values and relative inten- 
sities of various reflections are summarized in Table (3.1) and (3.2) for palladium 
and palladium silicide , respectively. Palladium film corresponds to f.c.c. phase 
with a=3.89A The diffraction data suggests that annealing of Pd/n-Si<lll> 
at 400-500"C invariably leads to silicide of composition PdaSi having FeaP struc- 
ture with hexagonal unit cell ; the parameters being a=13.055A , c=27.490A ^r]. 
Also , there is a tendency of preferred orientation , e.g. , at SOO^C , PdaSi grains 
assume preferred orientation such that their basal plane lies parallel to (111) of 
silicon substrate. (Notice appearance of 00.16 reflection and progressive increase in 
its intensity with annealing) . This means <111> Si is parallel to <00. 1> PdaSi. 
This result is consistent with the findings of Buckley and Moss [i'j. Further, at 
400°C and 450'’C , the intensities do not match with the ASTM data of poly- 
crystalline PdaSi and therefore indicate the presence of preferred orientation. The 
diffraction data suggests that (22.8) PdaSi lie parallel to (111) Si. With lapse of 
time , some readjustment occurs and grains assume other orientations as well and 
move towards true polycrystallinity ,i.e., PdaSi exhibits epitaxy in the begining but 
turns into polycrystalline disposition as it thickens. 

Palladium reacts with silicon in accordance with the equation 2Pd -f- Si — > 
PdaSi. Thickness d* of PdaSi formed can be estimated from 

j* Mpd^SiPPd 





4 ^ 



Figure 3.16: Diffractogram for the sample annealed at 200‘'C for differ- 
ent lengths of time. 
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Table 3.1: Diffraction peaks for Pd 


Pd 

Reflections 

hkl 






111 

200 

220 

311 

222 

as deposited 

29 

40.1 

47.1 

68.3 

82.5 

86.6 


dO) 

2.248 

1.930 

1.37 

1.171 

1.122 


R(I) 

85 

15 

30 

10 

8 

annealed 

29 

40.1 

" 47.1 

68.3 

82.5 

86.6 

at 200^^0 


2.248 

1.930 

1.37 

1.171 

1.122 

Ih 

R(I) 

80 

20 

40 

15 

5 


\29 

40.1 

47.1 

68.3 

82.5 

86.4 


dO) 

2.248 

1.930 

1.37 

1.171 

1.126 

lOh 

R(I) 

80 

18 

35 

10 

3 

annealed 

~29 

40.1 

47 

68.2 

82.4 

86.4 

at 400^0 

dai 

2.248 

1.933 

1.375 

1.173 

1.126 

0.5h 

R(I) 

180 

30 

60 

20 

10 


29 

40.1 

47 


82.4 

_ 


cKA) 

2.248 

1.933 

- 

1.173 

— 

6h 

R(I) 

80 

15 

- 

10 

- 


29 

40.1 

- 

- 

82.3 

_ 


d(A) 

2.248 

- 

- 

1.172 

— 

lOh 

R(I) 

30 


- 

5 

- 


29 

40.1 

- 

- 

- 

_ 


dC^) 

2.248 

- 

- 


_ 

16h 

R(I) 

20 

- 

- 

- 

- 

annealed 

29 

40.1 

- 

68.2 

82.5 

- 

at 450^C 

dCl) 

2.248 

- 

1.375 

1.171 

_ 

Ih 

R(I) 

50 

- 

28 

8 

- 


29 

40.2 

- 

68.3 

82.1 



d(^) 

2.246 

- 

1.373 

1.174 

- 

31i 

R(I) 

55 

- 

26 

7 

- 


29 

40.1 

- 

68.1 

82.5 



diX) 

2.248 

- 

1.377 

1.171 

_ 

5h 

R(I) 

40 

- 

20 

5 

- 


29 

40.2 

- 

- 

- 

- 


m 

2.246 


- 

- 

- 

7h 

R(I) 

20 

- 

- 

- 

- 

annealed 

29 

40.1 


68.5 

82.3 

- 

at SOO^C 

dO) 

2.248 


1.369 

1.172 

- 

0.5h 

R(I) 

50 

- 

30 

10 



29 

40.1 

- 

68.3 

82.5 

- 


dQd 

2.248 

- 

1.373 

1.171 

- 

1.5h 

R(I) 

44 

- 

28 

7 

- 


29 

40.2 

- 

68.8 

82.4 

- 


dlX) 

2.246 

- 

1.364 

1.1715 

- 

2.5h 

R(I) 

39 

- 

24 

5 

- 


29 

40.1 


68.2 

0 



dO.) 

2.248 

- 

1.375 


- 

5h 

R(I) 

26 

- 

20 

- 

- 
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Table 3.2: Difeaction peaks for Pd 2 Si 


PdaSi 

Reflections , bkl 



228 

408 

0016 

221j6 

2416 

6016 

as deposited 

29 

“ 

- 

- 


- 

- 



- 

- 

- 


- 



R(I) 

- 

- 

- 

• 

- 

- 

annealed 

26 


- 

- 

- 

- 

- 

at 200°C 


- 

- 

- 

- 

- 

- 

lOh 

R(I) 


- 

- 

- 

- 

- 

annealed 

29 

38.1 


- 

- 

- 

- 

at 400°C 

dih) 

2.362 

- 

- 

- 

- 


0.5h 

R(I) 

30 

- 

- 

- 

- 

- 


29 

38.1 

- 

53.7 

- 

- 

74.8 


dO) 

2.362 

- 

1.709 

- 

- 

1.270 

6h 

RW 

140 

- 

9 

- 

- 

24 


29 

38 

41.5 

53.8 

- 

70.3 

74.8 


d(h) 

2.366 

2.176 

1.703 

- 

1.339 

1.270 

lOh 

R(I) 

180 

32 

10 

- 

40 

26 


29 

38.1 

41.4 

53.5 

61 

70.4 

74.9 


dm 

2.362 

2.181 

1.713 

1.521 

1.337 

1.268 

16h 

R(I) 

200 

40 

12 

8 

50 

30 

annealed 

29 

38 

- 

53.3 

- 

70.3 

- 

at 450^C 

dm 

2.366 

- 

1.718 

- 

1.339 


Ih 

R(I) 

25 

- 

10 

- 

24 

- 


29 

38.1 

- 

53.3 

- 

70.3 

74.9 


dm 

2.362 

- 

1.718 

- 

1.339 

1.268 

31i 

R(I) 

35 

- 

12 


23 

5 


29 

38 

- 

53.4 

- 

70.4 

74.8 


dm 

2.366 


1.715 

- 

1.337 

1.270 

51i 

R(I) 

40 

- 

11 

- 

23 

10 


29 

38.1 

- 

53.3 

- 

70.3 

74.9 


dm 

2.362 

- 

1.718 

- 

1.339 

1.268 

7h 

R(I) 

60 


12 

- 

30 

16 

annealed 

29 


- 

53.7 

- 

70.3 

- 

at SOO^C 

dik) 

- 


1.706 

- 

1.339 


0.5h 

R(I) 

- 

- 

25 

- 

18 

- 


29 

38 


53.4 

- 

70.4 

- 


dm 

2.366 

- 

1.715 

- 

1.337 

- 

I. 5 I 1 

R(I) 

20 

- 

38 

- 

19 

- 


29 

38.1 


53.5 

60.9 

70.4 

- 


dm 

2.362 

- 

1.713 

1.521 

1.337 

- 

2.5h 

R(I) 

24 


42 

14 

19 

» 


29 

38.1 

41.5 

53.3 

60.9 

70.3 

74.9 


dm 

2.362 

2.176 

1.718 

1.521 

1.339 

1.268 

5h 

R(I) 

37 

5 

47 

15 

18 

10 




Results and Discussions 


47 


where Mp^^si and axe the molecular weight of Pd 2 Si and atomic weight of 
Pd , respectively, ppd, and pPd 2 Si their densities , dp^ thickness of Pd-film is deter- 
mined by optical interferometry as ~1050 ± 50A. Substituting the values of various 
paramaters as Mpd2S,=240.89 , Ap^=106.4 , Np^ = 2 , ppd = 12.02g/cm2 pPd^s, = 
9.462g/cm^ , dp^^g^ is estimated as 1510 ± 70 A. 

Wittmer and Tu have shown that the growth of Pd 2 Si on silicon wafer 
of <100> and <111> orientations is diffusion limited with activation energy (Eq) 
lying in the range 1.3-1. 5 eV. Also, Pd 2 Si follows parabolic growth law(i.e., d*^ = 
2(Dt)^/^, where d* is the silicide thickness, t is the annealing time and D(T)=Do 
exp (-Ea/kT) is the chemical interdiffusion constant; the value of Do is 0.48 or 
0.18 cm^/sec for pvdoped Si <100> or <111> orientation) and is independent of 
the substrate orientation.The rate limiting step is interdiffusion through the siUcide 
layer rather than the interfacial reaction . Further, the dominant diffusing species 
is palladium in Pd 2 Si via vacancies. 

These observations are at variance with the earlier reports of Cheung and 
Nicolet [51] who found diffusion constants to be in the range 10“^^- 10“^® cm^/s 
at 240"C and activation energy of Pd 2 Si formation as 0.9i0.1eV. Our value of 
interfacial diffusion constant lie in the range 2. 2-6. 3 x 10"^^ cm^/s at 400-500‘^C and 
eire closer to the later case. 

3.6 Reflectance measurements and SEM studies 

Fig. (3.20) , (3.21), (3.22) show reflectance(%R ) of Pd/n-Si system as such 
and after annealing for different lengths of time at 400, 450 and 500°C, respectively. 
Accordingly at 400°C, the reflectance spectra preserve the trend but decreases with 
time of piTinpflling . On careful examination, it was found that .patches of Pd 2 Si 
developed first at a ntunber of places , they enlarged to islands and eventually 
covered the entire area by coalescence. At 450 °C , the reflectance (%R) decreases 


Reflectance(%) 
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Wavelength(nm) 


Figure 3.20: Reflectance versus wavelength plots for the sample annealed at 400°C 
for different lengths of time 




Wavelength(nm) 


Figure 3.21: Reflectance versus wavelength plots for the sample annealed at 450°C 
for different lengths of time 


Reflectance(^) 



Wavelength(nm) 


Figxire 3.22: Reflectance versus wavelength plots for the sample annealed at 500°C 
for different lengths of time 
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initially , then increases upto 2h and eventually decreases at all wavelengths. In 
this case no patches , as observed at 400 °C , were found . Instead it is believed that 
lateral growth of Pd 2 Si occurs. While the first decrease can be attributed to initial 
growth of Pd 2 Si , the increase in the reflectance may be due to recrystallization and 
einc rgcnce of preferred orientations. As the time lapses , Pd 2 Si exhibits increased 
polycrystallinity (section 3.5) and so reflectance decreases as usual. At 500 “C , 
rc'flectance decreases with annealing time as basal plane (00.1) is predominantly 
lying parallel to the underlying (111) plane of silicon . But , after 1.5h , there 
is a slight increase in reflectance presumably due to other orientation developing. 
Eventually randomization occurs and we observe decrease of reflectance after 2.5h. 

The scanning electron micrographs are shown in fig(3.23) 
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Fimre 3 23' Scanning electron micrograpbsof Pd/n-Si Schottky contacts showi^ 
STstSnre (a)after annealing at 400»C for 12h (b) after anneabng at 450-C for 

41l (c) after annealing at 500»C tor 30 min. 


Chapter 4 
Simulation Studies 

4,1 I-V characteristics 

In order to understand the relative contributions of various current mechamsms, 
in Schottky diodes, I-V characteristics have been simulated for all of them using 
known i>arameters.For this, thermionic emission-diffusion, generation-recombination 
and tunneling processes are chosen.The thermionic emission-diffusion current through 
a Schottky diode under a forward bias of V volts is given by [4,7]. 

lie = (|^) - 1] (4.1) 

with saturation current Is given by 

I. = AiA'-T‘‘expf^) (4.2) 

where A** is the Richardson constant (1.12 x l{fATn~'^K~'^) and is the barrier 
height (for Pd^Si/n — Si contact, its value is ~ 0.735eV [11]) . A^ is diode area 
(for 1mm size dot A 4 = 7.85xl0“^m^) and k the Boltzmann constant is 1.38 x 
. With these parameters, Ig has been simulated using eqn(4.2) at various 
temperatures in the range 50K-300K. For each Is , in turn , thermionic emission- 
diffusion current {Ite) can be computed from eqn(l.lO) at various forward bias levels 
at the corresponding temperature. Figure 4.1 shows the results ln(I) vs (V) plots 
at various temperatures. Notice that saturation current Ig lies in the region ~ 
10 ~^j 4 — 10~®A at temperatures lOOK - 300K. 

The generation-recombination current through a Schottky diode under for- 
ward bias of V volt is given by [4]. 
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Figxire 4.1: Simulated la(I) vs V plots for the current due to thermionic emis- 
sion-diffusion process at different temperatures 
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, / 

= ( 4 - 3 ) 

with saturation current ciurrent 

Iio = qn,wAdl2r (4.4) 

wh('re, rjj is the intrin'jic carrier concentration, w is the width of the depletion region, 
Ad is the area of the junction(iuterface) andr is the carrier life time, n, is given by 

PI 


( 4 - 5 ) 

wh('re in* and m^ are the effective mass of the electron and hole respectively 
that have value 1.18m and 0.81m ; m being the mass of electron , is the band 
gap of silicon,k the Boltzmann constant and h is the Plank’s constant . The value 
of n, is determined at different temperatures with Eg = 1.12eV , w = 0.3/im, Ad = 
7.85xl0"^m^ and t = lO'^s [7]. Generation recombination current can be obtained 
from eqn(4.3) and eqn(4.4) at different temperatures (Fig.4.2). Notice that the 
saturation current is always 2-5 orders of magnitude less than the corresponding 
values obtained for the TED process (Fig.4.1). 

The tunneling current through the barrier under a forward bias of V volt is 

given by [4, 1X^(5] 


It = ItoexpiVlE') 


(4.6) 


with the saturation current Ito 
Ito ~ 

and every parameter 




k 


(4.7) 




(4.8) 
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Figure 4.2: Simulated In (I) vs V plots for current due to generation recombination 
process at different temperatures 
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E„ and E^o arc related such that 


Eo — Eoocoth 


qE(j^ 

~W 


(4.9) 


withi?^ given by 


47 r 771*65 


(4.10) 


where is the donar density and is permittivity of silicon. At donar con- 
centration Nfi ~ 1.35 X 10^^, m* = 1.18m and Cj = 11.9, E^ comes out to be 
1.88 X 10 eV. Using this value of Eoo,E' can be calculated at different temper- 
atures and in turn, ijo with A* = 1.12 x 10® 6;, = 0.763eU. Thus, tunneling 


current is obtained at different bias from eqn(4.6) for temperature range 50K-300K . 
Fig(4.3) shows variation of tunneling en rent with bias at various temperatures. The 
plots reveal very little variation of current with bias. Also,the current level drops 
down significantly with decrease in temperature. Fig (4.4) shows bias dependance 
of current more clearly at 300K. Obviously, tunneling current is always very low 
in comparison to current resulting due to TED and generation recombination pro- 
cesses.(See eg., Fig (4.5)). 

The ln(I)-V plots should be linear for all the three currents as per the govern- 
ing equations (4.1),(4.3),(4.6).It has however been observed that the experimental 
ln(I)-V plots are linear over several orders of current^higher biases.Such an effect is 
actually caused by the series resistance (Rs) present .In this situation , bias applied 
to the diode gets modified to (V-Ii?s)in equations(4.1)(4.3)(4.6) Saturation effects 
appear and the current can be estimated by iteration only as described in chapter 
3.Also,fitting of I-V data in TED current equation gives barrier height and j7.It 
is usually found that <j>ao decreases and r] increases with decrease in temperature 
and activation energy plot deviates from linearity at low temperatures.lt amounts to 
saying that the actual current is more than the contribution made by TED process 
itself. These anomolies are at first sight indicate a deviation from the thermionic 
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Figure 4.3: Simulated ln(I) vs V plots for the current due to tunneling process at 
different temperatures 
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emission-diffusion mechanism.The possible reason could be the simultaneous oper- 
ation of other processes viz., tunneling through the barrier or recombination in the 
dc'pletion region. 

The tunneling involves field emission and thermionic field emission through 
the barrier.Thc process applicable, however, depends upon the parameter 


Eoo = = 18.5 X (411) 

The field emission becomes important when kT /g,whereas thermionic 

field emission (TFE) dominates when E^o ~ kT/q and thermionic emission-diffusion 
if Eoo < kT/g [4]. In the present case Eoo = 0.188mV with Nn = 1.35xl0^^m"^ 

, rUg = 1.1 and = 11.9.This leads to qE^o/kT equal to 0.0071 at 300K and 0.027 at 
80K . Obviously , the condition prevailing in both the situations is Eoo ^ kTlq(^ 
0.259m 1^ and ~ 0.7mV at 300K and 80K, respectively) and hence the possibility of 
FE and TFE can be easily ruled out. The tunneling current as estimated from eqn. 
comes out to be very low ~ 10“^^ A at 300K and ~ 10~^^A at lOOK, whereas, the 
experimentally obtained values just at forward bias of O.IV are 10“®A and 10'^^ A 
at 300K and lOOK respectively. 

The second mechanism , i.e. , the generation and recombination requires (see 
,e.g., eq(4.3)) the ideality factor to be close to 2. But this value of t) is found at 80K. 
Below 80K,the q remains much smaller than 2.The saturation current calculated 
using eqn(4.4),comes out to be ~ 10“^^ Amp at 300K and 10“^° A at lOOK. Obviously 
the estimated current values for recombination process are excessively low compared 
to those actually observed in experiments.Chen et al [ 52 ] have pointed out that 
if recombination current exists, the barrier height evaluation on the basis of TED 
theory yields progressively lower values with decrease in temperature;for example, if 
Ir constitute 91% of the total current at (say) lOOK, barrier decrease expected is 
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0.02.This decrease is actually an order of magnitude lower than observed in Pd/n- 
Si Schottky diodesl^j*^. 

4.2 Temperature dependence of barrier height 

The second part of the simulation involves barrier height as a function of 
temperature. So far, barrier height is assumed to have a constant value at all 
temperatmes . However, the apparent barrier height of real Schottky contact shows 
an abnormal decrease as the temperature is lowered and discussed in length in 
chapter 3. It can be adequately accounted for if a spatial Gaussian distribution of 
barrier heights (with mean (pbo standard deviation <t ) is assumed to prevail at 
the interface. Accordingly , the apparent barrier height varies with temperature as 

^ ( 412 ) 

where and al are the mean barrier height and variance , respectively, at 
zero-bias. It may be pointed out that also corresponds to zero-bias only. In the 
above consideration (fbo and are taken to be independent of temperature. But, 
Werner and Guttler [3^] have shown both these parameters do vary linearly with 
temperature , such that 


<Pbo{T) = (f>bo{0) + OiiT 


(4.13) 


and 

<rl{T) = aliO)+a2T (4.14) 

where <fbo{0) and < 7 ^( 0 ) are the zero bias mean barrier height and variance , 
respectively , at absolute zero and m and aa their temperature coefficients . Now 
making substitutions into equation (4.12) , one gets, 
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= (4.15) 

Thus, contain terms amounting to direct and inverse temperature dependence. 
The aoove equation is used to determine (j}^^ in the temperature range 50K-350K 
for different values of and aa- For this data taken are (^bo{0K) = 0.8V and <r„(0) 
= 0.06V. The value of aj is varied from -0.3mV/K to +0.3mV/K and that of ag 
from -3(mV)VK to +3(mV)7K. 

Notice the coefficient Q 2 is for variance cr^. The corresponding values 
of temperature coefficient of standard deviation lie in the range -0.025mV/K to 
+0.25mV/K . Fig. (4.5) shows the variation of apparent barrier height as a function 
of temperature at various values oi when 02 is equal to zero. The situation with 
Oi = 0 corresponds to equation(4.12) itself. It depicts decrease of barrier height 
with decrease in temperature due to the presence of Gaussian distribution of barrier 
heights with mean(0i,c) and variance (al) . When <j)bo itself varies with temperature 
(i.e., 0 ;]= -3 to 3 mV/K and a%= 0 ) the extent of decrease is reduced or enhanced 
depending upon the sign of temperature coefficient (i.e., positive or negative). The 
difference is , however, more prominent at higher temperatures ,i.e., above 150K. 

The effect of 02 on the beirrier height is just opposite and shown in Figs. (4.6) 
and (4.7) for a]= -3mV/K and 3mV/K, respectively. Here , the extent of decrease 
is less or more if sign of 02 is negative or positive. Physically, it means that when 
Q 2 is negative (equivalent to narrowing of the distribution function, i.e., moving 
towards homogenization of barrier heights ) , the decrease of barrier height will 
be less than predicted by equation(4.12). On the other hand , when Q 2 is positive 
distribution function widens and amoimts to increase in the value of <To with 
rise in temperature . Consequently, there is effectively more reduction in ^bo then 
determined by equation(4.12). 
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Figure 4.6: Apparent barrier height versus temperattire plot with aj varying from 
-3 to 3 mV /K and zero 


Apparent barrier 
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Figure 4.7: Apparent barrier height versus temperature plot with mV/K and 
q 2 varying from -3 to 3 mV^/K 
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Figure 4.8; Apparent barrier height versus temperature plot with ais3 mV/K and 
a 2 varying from -3 to 3 mV^/K 



Chapter 5 
Conclusions 


• Pd/ii-Si(lll) Schottky barrier diodes can be prepared successfully by depositing 
a palladium layer of ~ 800-A thickness using an electron beam evaporation source 
in vacuum ~ 10~® mbar onto a thorougly cleaned silicon substrate. 

• characteristics of Pd/n-Si(lll) Schottky contacts can be explained on the 
basis of thermionic emission-diffusion (TED) mechanism . The ln(I) vs V plots 
over the temperature range (80K-300K) exhibit linearity over several orders of cur- 
rent and saturates due to effectiveness of voltage drop occuring across the series 
resistance (Rj) of the diode. 

• Several diode parameters , viz, barrier height , ideality factor , series resistance 
and saturation current can be determined by the fitting of the I-V data in the TED 
current equation using a computer program. While the zero bias barrier height 
decreases , ideality factor r) and series resistance increase with the decrease in 
temperature. Moreover , the changes in r) , and R* become quite significant at 
low temperatures. 

• The abnormal behaviour of Pd/n-Si Schottky diodes i.e. , decrease of zero 
bias barrier height and increase of ideality factor with decrease in temperature can 
be successfully explained on the basis of TED mechanism by incorporating the 
concept of barrier inhomogenities and describing them by a Gaussiein distribution 
of barrier heights. The values of parameters obtained are : zero bias mean barrier 
height (^6o) 0.93V and 0.63V , standard deviation (cto) 0.112 and 0.067V in the 
temperature range (165K-300K) and (80K-165K), respectively. 

• The <j)bo vs 1/T plot yields straight hne(s) to demonstrate the number of 
barrier distribution(s) , viz, single, double or multiple existing at the Schottky 
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contact giving values of both and ao . Similarly, inverse ideality factor (I/??) vs 
inverse t('mperature (1/T) plot displays linearity and it’s intercept at the ordinate 
and slope determine the bias coefficients (7 and of the mean barrier height and 
the standard deviation , respectively. 

• The annealing of Pd/n-Si <111> at 400-500°C invariably leads to silicide 
formation of the composition Pd2Si having a hexagonal structure (Fe2P type) with 
parameters a = 13.055A , c= 27.490A. Also, there is a tendency of preferred 
orientation in the begining but with lapse of time , some readjustment occurs and 
grains assume other orientations as well to exhibit true polycrystallinity. 

• Simulation studies suggest dominance of TED mechanism in the Pd/n-Si 
Schottky contacts. Also, negative temperature coefficients of mean barrier height 
{(fy) and variance (<t^) of the Gaassian distribution function (representing the bar- 
rier inhomogeneties at the Schottky contact ) are found to have opposite effects on 
the temperature dependence of apparent barrier height. While the former causes 
decrease in barrier height the later compensates the change by a fixed amount. 
The overall effect is that Gaussian distribution function shifts towards lower bar- 
rier height and at the same time becomes narrower with rise in temperature , i.e., 
effective reduction in barrier height occurs. 
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